Dissipative solitons and vortices in polariton superfluids 
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We examine spatial localisation and dynamical stability of Bose-Einstein condensates of exciton- 
polaritons in microcavities under the condition of off-resonant spatially inhomogeneous optical 
pumping both with and without a harmonic trapping potential. We employ the open-dissipative 
Gross-Pitaevskii model for describing an incoherently pumped polariton condensate coupled to an 
exciton reservoir, and reveal that spatial localisation of the steady-state condensate occurs due to 
the effective self-trapping created by the polariton flows supported by the spatially inhomogeneous 
pump, regardless of the presence of the external potential. A ground state of the polariton conden- 
sate with repulsive interactions between the quasiparticles represents a dynamically stable bright 
dissipative soliton. We also investigate the conditions for sustaining spatially localised structures 
with non-zero angular momentum in the form of single-charge vortices. 



I. INTRODUCTION 

Since the first observation of spontaneous exciton- 
polariton condensation in a microcavity significant 
effort in this vigorous field has been directed towards 
understanding the properties of superfiuid fiow in the 
novel non-equilibrium polariton superfiuid j^-iSj. Re- 
cently, considerable attention was drawn to the controlled 
creation of spatially localised collective excitations of the 
exciton-polariton Bose-Einstein condensate (BEC). Most 
notably, moving bright solitons have been predicted [6| 
and observed in a pioneering proof-of-principle experi- 
ment The properties of polariton solitons, superior 
to those of optical solitons in semiconductor cavity lasers, 
namely the picosecond response times and large nonlin- 
earities, suggest that polariton BEC can offer novel func- 
tionalities for information- processing devices [3]. The 
bright solitons observed in jSj] require a resonant excita- 
tion regime at nonzero in-plane momentum to make use 
of the dispersion properties of the lower-polariton branch 
that allow the polaritons to have a negative effective 
mass. In this regime, strong repulsive interactions be- 
tween the quasiparticles can be balanced by the effective 
dispersion to achieve localisation. Similar mechanism for 
spatial localisation via effective mass management is well 
explored for atomic BECs in optically induced band-gap 
structures [9]. In the absence of the spectral gap, how- 
ever, localisation by means of effective mass management 
is only possible in one spatial dimension. The question 
whether localisation of polariton BEC with repulsive in- 
teractions can occur in the regime of spontaneous con- 
densation at zero in-plane momentum (positive effective 
mass) and non-resonant excitation, remains open. 

In this paper we examine, both analytically and numer- 
ically, formation and dynamical stability of a steady-state 
BEC of exciton-polaritons in microcavities in the pres- 
ence of a spatially localised (Gaussian) optical pump. By 
emp loying the open-dissipative Gross-Pitaevskii model 
[10| describing an incoherently pumped BEC coupled to 
the exciton reservoir and successfully used in theoretical 



description of a number of significant experiments (see, 
G-g-, we analyse the mechanisms for creating and 

sustaining two-dimensional spatially localised structures, 
such as dissipative solitons and vortices. In addition to 
the trap- free case, we analyse the structure of the lo- 
calised states with the addition of a harmonic external 
potential that can be created, e.g., by engineered stress 
of the microcavity [ll|. In the latter case, the localisation 
of the repulsive polariton BEC is due to the harmonic 
confinement, as is the case with the analogous trapped 
atomic BEC with a positive scattering length. 

We show, analytically and numerically, that, even in 
the absence of external potentials localisation of the 
steady-state polariton BEC occurs due to the effective 
self-trapping created by the polariton fiows due to the 
spatially inhomogeneous off-resonant optical pumping. 
Spatial localisation of the condensate occurs due to the 
internal balance of superfiuid density fiows which is iden- 
tical to that responsible for supporting dissipative " anti- 
solitons" in optical systems. Such localisation of continu- 
ously self-defocusing solitons occurs despite the repulsive 
interactions between the quasiparticles, and is analogous 
to optical gain-guiding effect |12| . [Tsj in its reliance on 
continuous pumping. Thus a ground state of the polari- 
ton BEC represents a dissipative soliton which is spatially 
localised and dynamically stable, regardless of the pres- 
ence of the external potential. A harmonic trapping po- 
tential dramatically modifies the structure of the steady 
state due to the competition with the spatially inhomo- 
geneous pump (see Fig. [1]). 

A single vortex created in a localised steady state BEC 
by phase imprinting is similarly supported by the contin- 
uous inhomogeneous pump as a spatially localised dissi- 
pative vortex soliton. Akin to a vortex line in a trapped 
atomic BEC with a positive scattering length it is dy- 
namically unstable [14]. In the absence of potential, a 
vortex line spirals out of the polariton condensate and 
the condensate restores to its ground steady state. We 
show that the addition of fabricated (i.e. non-rotating) 
harmonic potential modifies the threshold of the optical 
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FIG. 1: Typical steady state density of a harmonically 
trapped polariton condensate in the regime of (a) narrow and 
(b) wide pump (see text). 



pumping required to sustain a steady state with an angu- 
lar momentum and leads to the possibility of a long-term 
survival of the vortex. 



II. MODEL 

We consider a spontaneously formed exciton-polariton 
BEC (fc = 0, lower-polariton branch) under the 
continuous-wave non-resonant excitation. The model 
first suggested in consists of a mean-field equation for 
the polariton condensate wavefunction and a rate equa- 
tion for the inhomogeneous density of the exciton reser- 
voir. We write it here in the form used in (^]: 



ih 



dt 
dt 



Vi + V{r, t) + i^iRuR - 7, 
2m z 



(1) 



where V{r,t) = C/(r) 0- Here \[> is the 
condensate wavefunction, hir is the exciton reservoir den- 
sity, P[r) is the inhomogeneous optical pump, and U{r) 
is the external potential. The critical parameters defining 
the condensate dynamics are the loss rates of the polari- 
tons 7c and reservoir excitons 7^, the stimulated scatter- 
ing rate i?, condensate coupling to the reservoir g^, and 
the coefficient gc quantifying the nonlinear interaction 
of polaritons. In what follows, we will consider a radi- 
ally symmetric external trapping potential U{r) — U^r'^, 
where t/o = in a trap-free case. 

The model can be re-written in dimensionless form by 
using the following characteristic scaling units of time, 
energy, and length: 



T = l/7c, E = h-f,, L : 



(2) 



where tolp is the lower-polariton effective mass. Here, 
we choose the parameters close to those of the experi- 



mental setup of Q, with ra^p — 10 
meVfj.m'^, gn = 2gc, 7c = 0.33 ps~ 



Ir 



= 6 X 10-3 
1.57c, and 



R = 0.01 fj,m? ps~^. In what follows we use the dimen- 
sionless (normalised) variables and parameters. 

We recall that, in the spatially homogeneous case, one 
can estimate the threshold pumping power at which the 
condensate appears in the system, as described in [To!]. 
This threshold is approximately determined by the values 
of parameters for which stimulated gain equals loss: 



th 
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R 



(3) 



The ratio of the polariton ric to exciton reservoir nji den- 
sities in the steady homogeneous state can then be de- 
fined as: 



nR 7fl, 



1), 



(4) 



where normalised pump is defined as P = P/Pth ^ 1- In 
general, even in a spatially inhomogeneous case, the typ- 
ical dynamical behaviour of the model Eqs. ([T]) displays 
rapid transition to the steady state of the polariton BEC 
for P > 1. 

Below we consider the properties of the steady state 
supported by an inhomogeneous (Gaussian) optical 
pump, P{r) — Pq exp(— r^/cr^), with or without an ad- 
ditional trapping potential. The spot size of the optical 
pump is considered to be small compared to the size of 
the sample, so that the local density approximation (l5j 
is not applicable. 



III. STEADY STATE OF THE CONDENSATE 
A. Theory 

For a general CW pump, the condensate wavefunction 
in the steady state can be looked for in the following gen- 
eral form: ^ = ^p{x,y) exp{—ifj,t), where fi is the chem- 
ical potential. If both the external potential and the 
pump are radially symmetric, we can re-write the model 
system in polar coordinates {x,y) — ?> {r,9) and reduce 
the problem to finding radially symmetric steady states 
ip{x,y) = '(/'(r) exp(jm^), where m is the phase wind- 
ing number [16] (topological charge of a vortex), and the 
ground state corresponds to m = 0. 

The steady state wavefunction ip{r) obeys the following 
stationary equation: 



Vl^j - 2V{r)'iP - i{Rn% 



7c)^ + = 0, (5) 



here = SV^r^ + (l/r)a/9r m^r^, V{r) = gc|V'P + 
gRTi'jj + U{r), and the steady state reservoir density is 
found from Eqs. ([T]) as: 



P(r) 



IR + R\i^\ 



1 ■ 



(6) 



All physical parameters in Eq. ([5]) are dimensionless, 
and under the adopted scaling ([2]), 7c = 1. However, we 
formally retain this parameter in subsequent formulas. 



3 



The stationary wave function can be further sep- 
arated into the real amphtude and phase: 'ipir) = 
$(r) exp[i<^(r)]. Following the analysis suggested in 
for the generalised complex Ginzburg-Landau equation, 
we can obtain asymptotic behaviour of amplitude and 
phase of the condensate wavefunction. In the absence of 
a trapping potential, U{r) = 0, the requirement of spatial 
localisation, P{r) -> 0, \tp{r)\ 0, imposed on both the 
pump and the condensate wave function at r ^ oo leads 
to the following asymptotic behaviour of the amplitude 
and the phase at large r: 

<i>(r) ^ Aexp{~p^r), (j){r) ^ B + p+r (7) 

1 /2 

where p± = [(/i^ + 7^/4)^/^ ± , and A,B are real 
constants. We note that, at large r, and irrespective of 
the vorticity to, the condensate decays slower than the 
pump field. 

Conversely, at r 0, the asymptotic behaviour is 
found to be: 

$(r) - Crl"l exp(-g_r2), 0(r) - q+r^ , (8) 

where g_ = (1/2)(1 + \m\)-^[V{0) - fi] and q+ = 
(1/4) (1 + \m\)~^[Rn%{0) - 7c], and for the ground state 
C = $(0) = $0- Note that for V{0) < the second 
derivative of the amplitude changes sign at r = 0, i.e. the 
condensate density displays a central "dip" in its spatial 
profile, as seen in Fig. [TJa). 

The presence of the harmonic trap U (r) — U'^r^ dra- 
matically modifies the asymptotic behaviour of the con- 
densate wavefunction for large r 18], so that: 

$(r) ^ AeM-Uy + (m - l/2)ln(2;7or)], (9) 
0(r) -S-|-7Jn(2i7nr). 

The asymptotic behaviour of a trapped polariton BEC 
at r is described by Eq. (|5]). 

We note that for a spatially localised (e.g., Gaussian) 
pump, P{r), the spatial localisation of the condensate 
in the absence of an added trapping potential is coun- 
terintuitive. Indeed, the effective potential V{r) formed 
by the repulsive nonlinearity (due to polariton scatter- 
ing) and interaction with the exciton reservoir is anti- 
trapping. Therefore the balance of nonlinearity and dis- 
persion responsible for nonlinear localisation in conserva- 
tive condensate systems cannot be achieved. Physically, 
the existence of the spatially localised steady state of 
the polariton condensate can be understood by exam- 
ining the internal flows of the dissipative polariton su- 
perfluid, analogous to the Pointing vector flows in dis- 
sipative optical systems [13, [ll]. Indeed, by defining 
the superfluid current density (flux) in the standard way 
j — Im(^'*V^), we can re-cast Eq. ([S]) in the form of 
coupled equations for the amplitude of the condensate 
wave function, $ and radial component of the stationary 
flux J = jr ^ n,cVr, whcrc Uc — ^^{r) is the condensate 
density, and Vr is the radial component of the flow ve- 
locity V = {d<j)/dr)er + {m/r)ee. The equations for these 



variables take the following form: 

i|;(rJ)-(i?n° -7c)$'-0, (10) 
- 2 [V{r) - ^] $ - K7$ = 0, 

where Vj = J^<i>~^ = {d(/)/dr)'^. The first equation of 
the system is the continuity equation for the stationary 
flow with source and sink. The first term in this equa- 
tion is simply the divergence of the fiux, D = V ■ j, which 
serves as a local measure of gain {D > 0) or loss {D < 0). 
The steady state exists if the generation of the polariton 
superfluid via continuous pumping is balanced by its dis- 
sipation. This condition can be formulated in terms of 
the flux as follows [T9| : 

/ D{r)rdr = 0. (11) 

Remarkably, this condition can be satisfied regardless of 
the sign of the nonlinear interactions in the systems, i.e. 
for both repulsive (as in the case of polaritons) and at- 
tractive nonlinearities. Nonlinear eigenstates for a very 
similar dynamical ID system with linear loss and spa- 
tially localised gain and without additional external po- 
tential were found in . The key feature of our system 
is similar: for the given parameters of gain and loss, the 
chemical potential, fi, is unique. 

The second equation in (ITU|) highlights the composition 
of the effective potential supporting the condensate wave- 
function as its eigenstate with the corresponding eigen- 
value /i. It is clear that the radial fiux that exists due 
to the spatially inhomogeneous phase of the condensate 
(j){r) ^ const, forms an attractive potential that can trap 
a spatially localised bound state even in the absence of 
the external trapping Uo — 0. In the linear limit (i.e. for 
very small condensates and / or nonlinearity) this effective 
trapping potential due to the flux is approximately given 
by: 




The effect of existence of the localised linear modes sup- 
ported purely by the localised gain is known as "gain 
guiding" in optics [l^ . 



B. Numerical results 

The existence of the localised steady state of the dis- 
sipative polariton BEC with the properties described in 
the previous section can be demonstrated by the numeri- 
cal simulation of the time-dependent model equation ([IJ . 
Fig. [2] shows the typical dynamical behaviour and tran- 
sition to the steady ground state (m = 0) of the polari- 
ton component for intermediate values of pumping power 
(Po = 3) in the trap-free regime. In addition, from the 
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numerical simulations we can extract and plot the follow- 
ing dynamical quantity: 



1 4 [V2|*|2 - 2|V*|2 - 4y|*|2] df 



(13) 



where S is the numerical integration domain. In the 
steady state limit this quantity coincides with the chemi- 
cal potential of the condensate Re[/i(t)] — >■ /i, Im[/i(i)] ^■ 
0. As can be seen from the dynamical simulations for 




Re(|i) 
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FIG. 2: Typical time evolution of (a) the condensate peak 
density and (b) the chemical potential for Po = 3, without a 
trapping potential. 

P > 1 (Fig. [2]), the steady state regime is reached quickly, 
and the chemical potential is well defined. The conden- 
sate wavefunction displays characteristic inhomogeneous 
amplitude and phase profile (/)(r) [seen in Fig. [3] with 
the limiting behaviour well described by Eqs (|7l8p . 

From Fig. it is seen that the inhomogeneous phase 
of the condensate wave function results in the nonzero 
flux of the polariton superfluid. Due to the spatially 
localised pump, the polariton superfluid is continuously 
generated in the core area of the steady state D > Q 
and dissipated on its wings D < Q. Thus the ground 
state of the strongly repulsive polariton BEC is a spa- 
tially localised dissipative soliton. Its spatial localisation 
is owed to the balance of the superfluid currents, and 
the internal structure of the currents is similar to that 
of a continuously "self-defocusing" dissipative "antisoli- 
ton" of the generalised complex Ginzburg-Landau equa- 
tion described in [l9|. 

The presence of a harmonic trap modifies the asymp- 
totic behaviour of the condensate and phase according 
to Eq. (|9]), as shown in Fig. [3l[b). Nevertheless, the 
structure of the internal flow within the ground state of 
the condensate remains qualitatively the same, as seen 
in Fig. [2Jd). In the trapped regime the internal flow 
of the polariton superfluid currents lead to significant 
distortions in the condensate density profile, nc(r). In 
particular, a wide pump with Pq ^ 1 leads to a well pro- 
nounced central peak [shown in Fig. [T] (b) for cr^ = 40 
and Po = 10] noted in numerical simulations of a sim- 
ilar model of polariton BEC [l^l and also hinted at in 
experimental observations of a trapped exciton-polariton 
condensate \Ya. 
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FIG. 3: Typical radial shape of the condensate ground state 
in (a) trap-free and (b) harmonically trapped case. Both 
amplitude, $(r), and phase, (l)(r), are shown in comparison 
to the Gaussian pump profile P{r)/Pth with the amplitude 
Po = P{0)/Pth ~ 6. Profiles of the radial flux J and its 
divergence D (scaled up by the factor of 10) for the ground 
state are shown for the (c) trap-free and (d) trapped cases, 
respectively. The harmonic potential of the depth Uq — 0.1 
is shown in (b) by the grey curve. 



IV. TRAPPED VS TRAP-FREE REGIME 

The ground m — state of the polariton condensate 
can be characterised by the dependence of chemical po- 
tential, fi on the parameters of the pump. In the absence 
of the trapping potential, this dependence is determined 
by the flux balance condition, and in the limit of low 
condensate densities the dependence of fi on the pump 
intensity is linear fi ^ Pq — P^/Pth- In the presence of 
the harmonic potential, this condition is modified by the 
trap, as seen in Fig. |4l 

For a fixed width of the pump, the cut-off value of 
Pq corresponds to the linear limit [■01 0, and devi- 
ates quite significantly from the value Pq = 1 deter- 
mined by the threshold behaviour ([3]) in the homoge- 
neous excitation case. The linear limit is determined by 
the cut-off for the ground state in the effective 2D po- 
tential formed by the combination of the external har- 
monic trap U{r), repulsive potential due to the spatially 
localised pump Vii{r) — gi^n'^, and the internal inhomo- 
geneous flux Vj(r) given by Eq. ([T2l) . This cut-off can 
be estimated from the condition that the effective poten- 
tial in the linear limit is trapping (attractive) rather than 
anti-trapping (repulsive). Thus, the bound state appears 
in the effective potential at the value of Pq given by the 
positive root of the equation: 



Pn' -2 



1 



4.97?, 1 



Po 



1 



4t^ 

v2 



= 0. 



(14) 



7cP a- J \ 7^ 

According to this formula, for the parameters in Fig. 
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FIG. 4: Dependence /^(Po) for the ground state m = without 
the harmonic potential (blue line) and in the presence of a 
harmonic trapping potential of the strength I7(r-) — U^r^. 
Vertical dashed lines indicate the cut-off value of Po for the 
ground state given by the linear limit ([14} . Marked points (a) 
and (b) correspond to the condensate profiles shown in Fig. 
El The width of the pump spot is = 30 in all cases. 
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FIG. 5: Dependence ^{(J^) for the ground state m = without 
and with a harmonic trapping potential (7(r) = (7o for Po = 
5. The vertical dotted line marks the cut-off value of the pump 
width for the trap-free ground state given by the linear limit 
expression (|14[l resolved with respect to . Marked points 
(a,b,c) correspond to the condensate profiles shown in Fig. (6] 



m the ground state appears in the trap-free system at 
Pq ~ 1-97. The presence of the harmonic trap with 
f/g = 0.1 lowers the threshold to Pq ~ 1.61. The pre- 
dicted tendency of the harmonic confinement to lower 
the threshold of the steady state formation compared to 
the trap-free case, agrees with the numerically calculated 
cut-off values in Fig. Above the threshold, the har- 
monic trap leads, for the same parameters of the pump, 
to the formation of the steady state BEC with stronger 
spatial localisation and higher peak density, as can be 
seen from comparison of Figs. Efa) and (b). 

For a fixed strength of the harmonic confinement, the 
width of the pump spot determines different trapping 
regimes. Intuitively, for an excitation spot much wider 
than the trap, the limiting behaviour should be described 
by a homogeneous pump approximation P{r) w Pq 
[isl [20j . In the opposite regime of a narrow excitation 
spot, the limiting behaviour should be close to that of a 
free condensate. Contrary to this intuitive picture, the 
dependence of the chemical potential /i on the width of 
the excitation spot, cr^, calculated numerically shows dra- 
matic differences for trapped and free condensates (Fig. 
[5]), and the condensate density profiles are strongly mod- 
ulated (see Fig. [6}. This is due to the presence of two 
competing spatial scales. One of them, ~ cr, is defined 
by the pump, and the other one ~ r^i?, is defined by the 
trapping potential and is given by the characteristic ra- 
dius of the wave function in the Thomas- Fermi limit. The 
latter is obtained by neglecting both phase and density 
gradients in Eqs. ([TU|. as the solution to the equation: 

rr2 2 , 1c9R f, ( '^\f\ /ii-n 

C^o'^Tf + — ^-Poexp ( — ^1=^. (15) 



In the regime of a wide pump, the condensate density 
can be reasonably well approximated by the radially- 
symmetric Thomas-Fermi profile with the radius t^f 
given by Eq.([T51). as can be seen in Fig. [Hla). In this 
regime, the ground state in the form of a dissipative soli- 
ton becomes dynamically unstable with respect to az- 
imuthal perturbations for Pq 3> 1, as established in [2(|. 
There it was shown that this instability can lead to rota- 
tional symmetry breaking whereby multiple vortex states 
enter the condensate. 

For a very narrow pump, <C t'j'p^ the Thomas- Fermi 
radius exactly coincides with the well-known expression 
for a conservative BEC with repulsive inter-particle inter- 
actions in a harmonic trapping potential: r^p — ii/Uq. 
In our numerical simulations we find that ttf ~ Tj-p for 
a wide range of values a < ttf ■ In this parameter regime 
the condensate density profile strongly deviates from the 
Thomas-Fermi profile [Figs. [6][b,c)], as neither phase nor 
density gradients can be neglected. The condensate ex- 
periences strong modulations of density with a notable 
density dip along the line D = separating the regions 
of loss and gain. The condensate peak density is also 
damped as the pump narrows down and the dip at r = 
appears for the values of a corresponding to the negative 
value of q- in We note that the point (b) in Fig. [5] 
corresponds to q_ = for Uq = 0.1. Under the condi- 
tions of a strong (Pq ^ 1) narrow pump the steady state 
is dynamically unstable and exhibits strong peak density 
oscillations. 
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FIG. 6: Radial shapes of the condensate ground states in the 
harmonically trapped case for Uq — 0.1, Po — 5 and the width 
of the pump spot corresponding to the points marked in 
Fig- m (a-c) Both amplitude, $(r), and phase, (pir), are 
shown in comparison to the condensate profile obtained in 
the Thomas- Fermi approximation, <1?tf (grey), (d-f) Profiles 
of the radial flux J and its divergence D (scaled up by the 
factor of 10) for the ground states shown in the panels (a-c). 



V. SINGLE VORTEX STATE 

Spatially localised dissipative vortices are similarly 
found as higher-order steady states of the system with 
?7i 7^ [m , both in a trap- free and harmonically trapped 
polariton BEC. Numerically, the system relaxes to a sin- 
gle steady state vortex when a phase factor exjp{im9) is 
imprinted onto a condensate "seed" at the initial stages 
of evolution. The dependence of the chemical potential 
/i on the pump amplitude Pq is plotted in Fig. [7] for a 
vortex with to = 1. As expected, the chemical potential 
and the cutoff for this excited state is higher than those 
of the ground state of the trap-free condensate. In con- 
trast, in a harmonic potential the vortex and the ground 
state become nearly degenerate (see Fig. [7]). 

The global phase of the vortex is imposed not only 
by the vorticity but also by the nonzero flux, and there- 
fore depends on both azimuthal (9) and radial (r) coordi- 
nates. Using the expression ([5]) for the radial dependence 
of the phase near the vortex core, we find that the lines 
of the constant phase follow a spiral trajectory described 
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FIG. 7: Top: Dependence /i(-Po) for the vortex state m = 1 
without and with the harmonic trapping potential U{r) = 
I/qT^ for = 30. Semi-transparent line indicates the corre- 
sponding dependence for a trapped m — state from Fig. |4] 
Bottom: (a) Condensate density and (b) phase for a dynam- 
ically unstable m — 1 stationary vortex state in a trap-free 
condensate at Pq = 8. (c) Metastable rotating vortex state 
supported by a harmonic trap (C/o=0.1) atPo = 6 



by the equations: 



Xs (9) - y/9/^ cos 9, y, (9) - y/9/^ sin 9 

This spiral phase structure is shown in Fig. ^h) for the 
numerically found localised vortex state in Fig. [71[a), 
without the trapping potential. It is a signature of the 
polariton BEC vortex observed in the experiments |2j, 
and also appears as a characteristic feature of vortex 
states in other dissipative systems [13, [H, [1^ . 

The ground state (to = 0) of the polariton condensate 
is dynamically stable in the absence of a trap. In con- 
trast, we find that a single vortex created in a localised 
steady state BEC by phase imprinting is dynamically un- 
stable and, in the absence of potential, spirals out of the 
condensate, as predicted for the case of trapped dissipa- 
tive BEC of alkali atoms 

MM- As can be seen from 

comparison of Figs. |4] and [71 the addition of a station- 
ary (i.e. non-rotating) harmonic potential dramatically 
modifies the threshold of the optical pumping required to 
sustain a steady state of a polariton BEC with an angu- 
lar momentum and leads to the possibility of a long-term 
survival of the vortex in the form of an eccentric rotat- 
ing state shown in Fig. [TJc). The detailed study of the 
stabilising effect of the trapping potential on the vortex 
dynamics and the dynamics near the onset of rotational 
symmetry breaking [20] is beyond the scope of this study. 
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However, we have confirmed that a wide pump (corre- 
sponding to small density gradients in the harmonic trap) 
tends to stabilise the vortex state. 

VI. CONCLUSIONS 

We have characterised the 2D stationary regime of the 
polariton BEC in the framework of an open-dissipative 
Gross-Pitaevskii model coupled to an exciton reservoir 
and described, both analytically and numerically, proper- 
ties of the ground and excited states (vortices) depending 
on the pump and trapping parameters. We have shown 
that the spatial localisation of the condensate, even in 
the absence of a trapping potential, is supported by the 
balance of the internal superfluid flows established by an 
inhomogeneous nonlinear gain due to optical pump and 
linear loss due to the decay of the polaritons. The ground 
state of the condensate can therefore be described as a 
continuously self-defocusing dissipative soliton. We have 
also investigated localised dissipative vortex states and 
have shown that they can display metastable behaviour 
in non-rotating trapping potentials. 



Finally, we note that the open-dissipative model with 
the inhomogeneous pump used here and successfully em- 
ployed for theoretical description of experiments with mi- 
crocavity polaritons, has many similarities to optical sys- 
tems with localised gain landscapes [13,^]. These and 
similar studies suggest that the 2D stationary regime 
of polariton BEC can potentially display a rich variety 
of localised states which are not exhausted by radially 
symmetric configurations. In particular, the experiments 
with polariton BEC excited by elliptical optical pump 
have demonstrated pattern formation due to the appear- 
ance of standing waves (25| . Furthermore, the recent 
experiment with polariton BECs created by two pump 
spots [2^ demonstrates novel possibilities arising from 
interaction of two dissipative solitons with variable sepa- 
ration. The clear understanding of the structure of phase 
gradients underpinning the existence of steady-state po- 
lariton condensates presented in our work enables us to 
gain an immediate insight into the properties and out- 
comes of such interactions. The possibility to create 
a rich variety of localised states by varying the spatial 
properties of the off-resonant optical pump are currently 
under investigation and will be reported elsewhere. 
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